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D ' We investigate the distribution of the spacings of adjacent eigenvalues of the lattice 

Dirac operator. At zero chemical potential fi, the nearest-neighbor spacing distribution P(s) 
• follows the Wigner surmise of random matrix theory both in the confinement and in the 

deconfinement phase. This is indicative of quantum chaos. At nonzero chemical potential, 
the eigenvalues of the Dirac operator become complex and we discuss how P(s) can be 
I . defined in the complex plane. Numerical results from an SU(2) simulation with staggered 

fermions in fundamental and adjoint representations are compared with predictions from 
lyv non-hermitian random matrix theory, and agreement with the Ginibre ensemble is found for 

y—{ i M~0.5. 
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q ' §1. Introduction 

^ : 

The spectrum of the QCD Dirac operator, both in the continuum and on the 
lattice, has several universal properties. By "universal" one means "independent of 
the details of the dynamics", e.g., independent of the precise values of the simula- 
tion parameters on the lattice. Such universal features can be described by random 
matrix theory (RMT). In this contribution, we are concerned with the eigenvalue 
fluctuations in the bulk of the spectrum on the scale of the local mean level spacing, 
measured by the distribution P(s) of the spacings s of adjacent eigenvalues. We will 
consider gauge group SU(2) and staggered fermions in the fundamental representa- 
tion which are related to the chiral symplectic ensemble of RMT. At zero chemical 
potential /x, all Dirac eigenvalues are purely imaginary, and it has been shown in 
lattice simulations that P(s) agrees with the Wigner surmise of RMT, 

D , , 262144 4 _6± s 2 

both in the confinement and in the deconfinement phase. ^ This result implies that 
the Dirac eigenvalues are strongly correlated, and is indicative of quantum chaos, 
according to the conjecture by Bohigas, Giannoni, and Schmit. 2 ) In contrast, quan- 
tum systems with uncorrelated eigenvalues, corresponding to classically integrable 
systems, obey a Poisson distribution, 

Pp(s) = e~ s . (1-2) 

Further, we present a first analysis for gauge group SU(2) with staggered fermions 
in the adjoint representation being related to the chiral orthogonal ensemble with 
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the Wigner surmise 



p w{s) = - se 4 s 



(1-3) 



While fermions transforming with respect to the fundamental representation de- 
scribe quarks, fermions transforming according to the adjoint representation can be 
interpreted as gluinos. Thus the latter case with 



1 



(1-4) 



constitutes a lattice version of supersymmetric Yang-Mills theory. 3 ) 

We focus on nonzero chemical potential, ^ / 0, where the lattice Dirac matrix 
generalizes to 



M XiV (U,iS) = — ^2 [Uv(x)rjv(x)6 y>a> {i, - h.c] 



u=x,y,z 



+ 



2a I 



y,x-t 



(1-5) 



with the Kawamoto-Smit phases n. The eigenvalues of this matrix are complex. 
In this case, P(s) represents the spacing distribution of nearest neighbors in the 
complex plane. For each eigenvalue zq one has to find the eigenvalue z\ for which 
s = \z\ — zo\ is smallest, with a subsequent average over zq. This definition assumes 
that the spectral density is constant over a bounded region in the complex plane 
(and zero outside). Generally, this is not the case so that an unfolding procedure 
must be applied; see Sec. El However, the spectral density of the so-called Ginibre 
ensemble of RMT, where real and imaginary parts of the eigenvalues have the same 
average size, is constant inside a circle and zero outside, respectively. 4 ) In this case, 
P(s) is given by 5 ) 

P G (s) = cp(cs) , (1-6) 



with 



n 



p(s) 
x /m\ 



2 s lim 

N— »oo 



N-l 



Y[ e n (s 2 ) e~ 



N-l 



■In 



^ nle n (s 2 ) 



m=o * and c = J °° ds sp(s) = 1.1429... 

In contrast, the Poisson distribution in the complex plane, representing uncorrelated 
eigenvalues, becomes 

Pp{s) = \se-^ 2 . (1-7) 

This should not be confused with the Wigner distributions (|1T|) and (jl-H[) . In the 
following, we will study the Dirac spectrum on the lattice at various values of /j, ^ 
and compare the resulting P(s) with Eqs. ()1-6|) and (|1-7|) . 



§2. Analysis of complex spectra 



We report on simulations done with gauge group SU(2) in the confinement re- 
gion; specifically we studied a 6 4 system at = 1.3 with 2 flavors of staggered fermion 
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Fig. 1. Complex eigenvalues of the Dirac operator at various values of p for a typical equilib- 
rium configuration of two-color QCD with staggered fermions in the fundamental and adjoint 
representation (horizontal axes = real parts, vertical axes = imaginary parts, in units of 1/a). 
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with bare mass ma = 0.07 in the fundamental representation, 6 ) and a 4 4 system at 
j3 = 2.0 with 1 flavor of staggered fermion with ma = 0.1 in the adjoint represen- 
tation (corresponding respectively to Nf = 8 and 4 continuum flavors). For these 
systems the fermion determinant is real and lattice simulations become feasible. 7 ) 
We sampled a few 100 - 1000 independent configurations for each p-vslue. 

The eigenvalue spectrum is shown in Fig. ^ for three different values of p. The 
size of the real parts of the eigenvalues grows with p as expected. While the spec- 
trum is confined to a bounded region in the complex plane, the spectral density is 
certainly not constant. Therefore, the spectrum has to be unfolded. In one dimen- 
sion, unfolding is a local rescaling of the eigenvalue density such that the density 
on the unfolded scale is equal to unity. To unfold spectra in the entire complex 
plane we proceed as follows. The spectral density has an average and a fluctuating 
part, p(x,y) = p sv (x,y) + pa(x,y). Unfolding means to find a map z' = x' + iy' = 
u(x, y) + iv(x, y) such that p scv (x', y') = 1. Since the probability has to be invariant, 
Pa,v( x ' ,y')dx'dy' = dx'dy' = p av (x,y)dxdy, we find that p av (x,y) is the Jacobian of 
the transformation from (x,y) to (x',y'), p av (x,y) = \d x ud y v — d y ud x v\ . Choosing 
y' = v(x,y) = y yields p av (x,y) = \d x u\ and, thus, x' = u(x,y) = f^dtpgy&y) . 
Essentially, this is a one-dimensional unfolding in strips parallel to the real axis. 8 ) 
For a fixed bin in y, p av (x,y) is determined by fitting p(x,y) to a low-order polyno- 
mial. P(s) is then constructed from the constant unfolded density as explained in 
Sec. n normalized such that ds s P(s) = 1. 

Several remarks are in order, (i) We have checked spectral ergodicity. If only 
parts of the spectral support are considered, the results for P(s) do not change, (ii) 
If the spectral density has "holes" , we split the spectral support into several convex 
pieces and unfold them separately. This is justified by spectral ergodicity. (iii) 
Unfolding each spectrum separately and ensemble unfolding yield the same results 
for P(s). (iv) The results for P(s) are stable under variations of the degree of the 
fit polynomial and of the bin sizes in x and y. 

§3. Results and discussion 

Our results for P(s) are presented in Fig. |2| As a function of p, we expect to 
find a transition from Wigner to Ginibre behavior in P(s). This was clearly seen in 
color-SU(3) with Nf = 3 flavors and quenched chemical potential, 9 ) where differences 
between both curves are more pronounced. For the symplectic ensemble of color- 
SU(2) with staggered fermions in the fundamental representation, the Wigner and 
Ginibre distributions are very close to each other and thus harder to distinguish, but 
they are reproduced for p = and p = 0.4, respectively. For the orthogonal ensemble 
with staggered fermions in the adjoint representation this transition is more drastic. 

Increasing p > 1.0, the lattice results for P(s) deviate substantially from the 
Ginibre distribution and can be interpreted as Poisson behavior, corresponding to 
uncorrelated eigenvalues. (In the Hermitian case at nonzero temperature, lattice 
simulations only show a transition to Poisson behavior for (3 — > oo when the physical 
box size shrinks and the theory becomes free. 1 )) A plausible explanation of the 
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Fig. 2. Nearest-neighbor spacing distribution of the Dirac matrix with fundamental and adjoint 
fermions for various values of /i. The histograms represent the lattice data. The curves corre- 
sponding to the Wigner, Ginibre and Poisson distribution are inserted. 
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transition to Poisson behavior is provided by the following two (related) observations. 
First, for large \i the terms containing in Ea. ()l-5j) dominate the Dirac matrix, 
giving rise to uncorrelated eigenvalues. Second, for large [i the fermion density on 
the finite lattice reaches saturation due to the limited box size and the Pauli exclusion 
principle. 

In conclusion, we have investigated the nearest-neighbor spacing distribution of 
the lattice Dirac operator for two-color QCD with fermions in the fundamental and 
adjoint representation switching on a finite chemical potential. We find a transition 
from a Wigner to a Ginibre distribution around the phase transition for both types 
of fermions. This means that quantum chaos persists deep into the high-density 
phase. 
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